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^ ■ Abstract 

i 1 The vast majority of works on option pricing operate on the assumption of risk 

neutral valuation, and consequently focus on the expected value of option returns, 
and do not consider risk parameters, such as variance. We show that it is possible 
to give explicit formulae for the variance of European option returns (vanilla calls 

c +h | and puts, as well as barrier options), and that for American options the variance can 

be computed using a PDE approach, involving a modified Black-Scholes PDE. We 
show how the need to consider risk parameters, such as the variance, and also the 
l | probability of expiring worthless (PEW), arises naturally for individual investors in 

' options. Furthermore, we show that a volatility smile arises in a simple model of 

ly-) ■ risk-seeking option pricing. 

-xt- : 

' 1 Introduction 

o ■ 

7-h ' Options are utilized by different investors for different goals. On the one hand options are 
often used for hedging, eliminating the risk of extreme downside events at the expense of 
forfeiting (some of) the gain of upside events. On the other hand, options are used for 
speculation, to increase the leverage of an investor that has strong convictions about the 
direction of the market and is willing to expose herself to substantial risk, with the potential 
of very large gains if the scenarios she anticipates are realized. Other uses are discussed in 

Regardless of the strategy in which they are being used, an option should be considered 
as a risky asset, with a distribution of returns determined by the distribution of returns of 
the underlying asset, but typically far more extreme, usually with significant probability of 
total loss or a gain of hundreds of percent. In this context it is bizarre that the most common 
piece of quantitative data given about options is the Black-Scholes price, which represents 
the current value of the expected value of the return of the option (under suitable assumptions 
about the distribution of returns of the underlying asset). Given that the distribution of 
the returns of the option is extreme, knowing its expected value is of limited significance, 
except maybe for the largest investors who are making many purchases and sales, and 
even these investors should be interested in knowing the relevant law of large numbers for 
their situation, and must take into account strong correlations. More information about 
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the distribution of option returns is surely needed in order to make rational investment 
decisions. 

Rather remarkably, we have not succeded in finding any literature examining the full 
distribution of the returns from options, even on the assumption of a lognormal distribution 
of returns of the underlying asset, beyond acknowledgment that this distribution is extreme, 
and, of course, calculation of its expected value. In part this can be attributed to the fact 
that seminal papers such as Black and Scholes [TJ, Merton [T7], Cox, Ross and Rubinstein 
[6], all emphasize that any deviation of an option price from its risk-neutral valuation creates 
an arbitrage opportunity that cannot exist in a rational market. In addition, much literature 
focuses on the variation of returns from options due to lack of complete knowledge about the 
distribution of underlying assets, their rates of return and (particularly) their volatilities, 
see for example Coval and Shumway [5], Liu and Pan [16] and Goyal and Saretto [12]. 
It is not the intention of this paper to question the rationality assumptions behind risk- 
neutral option valuation. The thrust of this paper is to show that mathematically it is 
possible to compute different parameters describing the distribution of returns from an 
option, assuming a lognormal distribution of returns of the underlying asset, with known 
parameters. We focus on two (related) measures of risk: the variance and the probability 
of expiring worthless (PEW). In the case of vanilla European call and put options these 
parameters are straightforward to compute analytically, and we do this in section 2. Analytic 
expressions can also be given for European options with barriers and we illustrate this in 
section 3. For American options, things are more complicated. Analytic expressions cannot 
be found, but the moments of the distribution from which the variance and PEW can be 
derived satisfy modified Black-Scholes partial differential equations (PDEs) giving numerical 
schemes for computation. This is described in section 4. 

In section 5 we give some simple examples of the use of the risk parameters we compute 
in individual investment decisions. In section 6, we use the variance parameter to build a 
simple model describing risk-averse or risk-seeking market pricing of options; we show that 
such a model is equivalent to the standard model with a volatility smile. Finally, section 7 
contains closing comments and suggestions for further research. 

2 European options 

Assuming an underlying with a lognormal distribution with parameters r, a and current 
price 5*o, the current value of the return on European call and put options with strike K at 
time T are 



where in both formulas Z is a standard normal variable. A simple computation gives the 
expected values 
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These formulae are the Black-Scholes prices for calls and puts respectively. Note both 
u = E[C] and u = E[P] satisfy the Black-Scholes PDE 

du 1 2 2 d 2 u du 

dr = r s °dsi +rSo ds- - ru - (1) 

It is a straightforward calculation to find the variance of C and P. We find 

E[C 2 } = Sy 2T <S>(-d + 2aVT) - 2KS e~ rT <b{-d + aVT) + K 2 e~ 2rT ^{-d) , 
E[P 2 ] = K 2 e- 2rT $(d) - 2KS e- rT $(d - aVT) + S V 2T $(d - 2a VT) . 

The variances are then computed using Var(C) = E[C 2 ] — (E[C]) 2 and Var(P) = E[P 2 ] — 
(E[P]) 2 . By the Feynmann-Kac formula, both v = E[C 2 ] and v = E[P 2 ] satisfy the modified 
Black-Scholes PDE 

dv 1 2 2 d 2 v dv 

&r = 2 aS °ds 2+rSo ds- - 2rv - (2) 

The factor 2 before the last term arises as the discount factor in C 2 and P 2 is e~ 2rT . The 
relevant final condition and boundary conditions for the modified Black-Scholes PDE (for 
both call and put) are obtained be requiring that the variance should vanish as the final 
time is approached, and as S — > 0, oo. 

In greater generality it is also possible to compute formulae for the nth moments E[C n ] 
and E[P n ], n = 3, 4, The nth moments satisfy the modified Black-Scholes PDE 

dv 1 2 r-2 d 2v n 9v 

dT = 2 aS °dSl + rSo dS- - nrV 

due to the discount factor e~ nrT . There does not appear to be an explicit closed form for 
the moment generating functions of C and P, but it can be shown that these also satisfy a 
suitably modified Black-Scholes PDE. 

Both the random variables C and P have a mixed discrete-continuous distribution, with 
a finite probability of being zero, and a continuous range of values above zero. For such 
distributions there are two possible ways of defining the zeroth order moment, either as 1 or 
as the probability to be nonzero. For options, the probability of expiring worthless (PEW) 
is a quantity of interest. Clearly we have 

PEW(C) = , 
PEW(P) = . 

Thanks to the relation of these quantities with zeroth order moments, they too solve a 
modified Black-Scholes PDE 

9w 1 2Q 2 d2w . c dw 

dT = 2 (lS °dS 2+rSo dS- - (3) 

In Figure [1] we show the variance and PEW of a put option with r = 0.1, a = 0.15, 
K — 1 as a function of Sq and for a variety of different times T. As expected, the variance 
is small for both low and high values of Sq, and in general increases with T. The location 
of the maximum variance falls as T increases. The PEW tends to a step function for small 
times. 
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Figure 1: Variance and PEW for a European put option with r = 0.1, a = 0.15, K = 1 



3 A down-and-out put option 

The variance and the PEW can also be computed analytcially for European barrier options. 
We illustrate this in this section by considering the case of a down-and-out put; that is, a 
European put option with the usual 5 parameters r, a, T, K, So and an additional "barrier" 
B < mm(K, So); if at any time during the life of the option the price of the underlying 
falls below B, then the value of the option is knocked to zero (irrespective of how the price 
develops after hitting the barrier). 

We use the following result proven in pTj : if Y(t) —W(t) + bt is a Weiner process with 
drift then 



P(min < t < T F(t) < a,x < Y(T) < x + dx) = y^ ex P [ 2 T j d% (4) 



where 

1 a > x or a > 



A 



exp ( 2a{x T a) ) otherwise 



Note that if A = 1 the RHS in fl4]) is simply P(x < Y(T) < x + dx). If a > or a > x 
then the requirement that min <t<TV(t) < a is automatically satisfied, as Y{t) starts from 
and ends at x. This is why A = 1 in these cases. Remarkably the correction factor A does 
not depend on b. It is straightforward to verify that the probability density in 01]), in both 
regions, satisfies the forward Kolmogorov equation 



The current value of a put option can be written in the form 

P = e- rT (K-S e° Y ^) + 

where Y(T) is a Weiner process with drift b — - — |. To avoid hitting the barrier we need 
mm < t < T Y(t) > a where a = ^ In (-§^J < 0. To obtain a nonzero return from the option 
we need Y(T) < ^ In Putting this all together we deduce that in the presence of the 



dx 



knock-out barrier at B 

PEW{P) . 

where we have written 

a = — In I — ) and b = . 

a \S J a 2 

Throughout we are assuming B < mm(K, So)- All the integrals are Gaussian, and the 
calculations are straightforward though tedious. The final results are as follows: 

PEW(P) = 1 - U(Q ) - S(Q0 + ^ 1 ($(Q 2 ) - $(Q 3 )) J , (5) 



where 



0l = 1 ^ftUf.-i^T 



5, 



+i 

-5o|$(g4)-$(Q 5 )+(^)^ (HQe) - HQr)) I (6) 



where 
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E[P 2 ] = KV 2rr |$(Qo)-W+(^r $(Q 3 )) 

-2/v,s;,-^ I $(Q 4 ) - $(g 5 ) + (ci,(0 (;) - ci> ( () 7 )) 



5', 



where 
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+Sy 2T ( $(Q 8 ) - $(Q 9 ) + (-) ^ ($(Q 10 ) - $(Qh)) 1 (7) 
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Formulae ([5]), ([6]) and ([7]) are of limited importance in themselves. Their significance lies 
more in the fact that, as can be readilly verified using a symbolic manipulator, Formulae 
(TJP, and |7|) satisfy the (modified) Black-Scholes equations (T^j, (QJ) and (d)) respectively. 
As we shall see in the next section, this can provide a numerical procedure to determine the 
quantities we need even when analytic expressions are not available. 

We conclude this section by presenting plots of the PEW, expected value and variance 
of a down-and-out put option with the same parameters as the European option considered 
in section 2, but with a barrier at B = 0.5. See Figure [21 For Sq well above the barrier the 
results are similar to those for the option without a barrier. But as So decreases towards the 
barrier, the PEW increases towards 1, the expected value drops and the variance increases 
rapidly, before dropping again to when So = B. 



4 American options 

For the case of American options there is no analytic expression for the Black-Scholes price, 
and we do not expect there to be analytic expressions for the variance or the PEW. (We 
limit the discussion here to American puts on non-dividend paying stocks; if no dividends 
are paid then for call options there is no added value in the early exercise feature [2TJ. The 
discussion can be extended to the case of continuous-dividend paying stocks.) The preferred 
approach for pricing American options (on a single underlying asset) is by numerical solution 
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PEW Expected Value Variance 




Figure 2: PEW, expected value and variance for a down-and-out barrier put option with 
r = 0.1, a = 0.15, K= 1, B = 0.5. 



of a free boundary problem for the Black-Scholes PDE. For an option with strike K and 
expiration T on an underlying with volatility a, the Black-Scholes price u(s, t) (as a function 
of the underlying price s and time t) should satisfy the PDE 

du 1 9 9 d 2 u du , , 

on the domain < t < T, s > b(t), where b(t) is the unknown early exercise price. u(s,t) 
should satisfy the Dirichlet boundary conditions 

lim u(s, t) = , 

s— >oc 

u(b(t),t) = K-b(t) 

and the final time condition 

u(s,T)=0, s>b{T). 
The early exercise price b(t) is determined by the smooth pasting condition 

u s (b(t),t) = -1 

along with the final time condition b(T) = K. See [21] for a full exposition. 

Completely analogously to our findings in the previous 2 sections, the expectation v(s, t) 
of the square of the current value of the option should obey the modified Black-Scholes 
equation 

dv 1 2 2 d 2 v dv n , n . 
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on the same domain, with boundary conditions 

lim v(s,t) = , 

s— >oo 

v(b(t),t) = (K-b(t)f 

and final time condition 

v(s,T) = 0, s>b(T). 

These conditions enforce vanishing variance at all boundaries. Note that there is no reason 
to expect the variance to have zero derivative at the early exercise boundary — thus it need 
not paste smoothly onto the variance below the early exercise boundary, which vanishes. 
Similarly, the PEW w(s,t) should obey 

dw 1 o 9 d 2 w dw , . 



with boundary conditions 



lim wis, t) = 1 

s—>oo 

w(b(t),t) = 



and final time condition 



w(s,T) = l. s>b(T). 

Note that for the PEW there is a discontinuity in the boundary data at t = T, s — B(T). 

Note that the early exercise price b(t) is determined by the free boundary problem for 
the standard Black-Scholes equation (JHJ), and this is then taken as input to solve the relevant 
problems for iQ and (fTUl) . In practice we solve all 3 PDEs using a front-fixing technique 
, [H]. Define new dimensionless coordinates y, t via 



and rescale the unknown functions u, v, w, b via 



u(s, t) = KU{y, t) , v(s, t) = K 2 V{y, r) , w(s, t) = W{y, r) , b{t) = KB{r) 
The system of equations becomes 



dU 


d 2 U 


!h 


dy 2 


dV 


d 2 V 


a7 


dy 2 


dW 


d 2 W 


dr 


dy 2 



all on the fixed domain < r < r max = \o 2 T , y > 0, with initial conditions 

U(y, 0) = V(y, 0) = , W(y, 0) = B(0) = 1 , 
Dirichlet boundary conditons 

U(0, r) = 1 - B(t) , V(0, t) = (1 - B(t)) 2 , W(0, t) = 
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Figure 3: PEW, expected value and variance for an American put option with r = 0.1, 
a = 0.15, K — 1. Including a blow-up of the intersection region for the PEW. The PEW 
and variance go to zero on the early exercise boundary. 



and 

lim U(y, t) = lim V(y, r) = , lim W (w, r) = 1 , 
and finally the Neumann condition 

^(0,r) = - J B(r) . 

We solve this system using the Crank-Nicolson scheme to advance U, V, W; the simplest first 
order discretization is used to approximate B'(t), and a 3-point, second order approximation 
is used for U y (0, r). At each time step the Newton- Raphson method is used to update B in 
a manner that the Neumann boundary condition holds. Results are displayed in Figure E] 
for an option r = 0.1, a = 0.15, K = 1. By construction the variance and PEW go to zero 
on the early exercise boundary. The variance apparently increases uniformly with T (for 
fixed s), which makes good sense — but note this was not exactly the case for the European 
option (see Figured]). As T tends to zero the PEW tends to a step function, as expected. 

A few notes on numerical aspects of the above calculations: The discontinuity in the 
derivative of the PEW at the early exercise boundary did not interfere with convergence of 
the Crank-Nicolson method, but the usual 0(h 2 ) error behavior was reduced to 0(h). The 
variance calculation was verified using a Monte-Carlo technique; once the early exercise curve 
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has been calculated, standard Monte-Carlo techniques can be used to get the distribution 
(and in particular the mean and variance) of the return from the option. These results 
confirmed the PDE calculations, but a "barrier correction" following jl] was needed to 
compensate for discrete-time observation of whether the stock hits the early exercise price 
in the Monte-Carlo simulation. There are also Monte-Carlo techniques for pricing American 
options that do not depend on knowing the early exercise price, for example the method of 
[21 [3] uses a tree approach to compute two estimates which are lower and upper bounds for 
the option price. The variances of the lower and upper bounds computed in this approach 
were found to be substantially lower than the option variance computed by PDE methods. 
This is perfectly reasonable, the choices made to obtain bounds eliminate a lot of the 
variation of the option price. 

5 Variance and PEW in Investment Decisions. 

Consider an investor buying a put option on a stock. The investor faces the problem of 
selecting from a universe of options indexed by their strike K and time to expiration T, 
with prices increasing as functions of both K and T. Assuming the options are priced by 
the classical Black-Scholes formula (for European puts), how is the investor to make this 
decision? 

The expected rate of return does not depend which options the investor buys (and for 
that matter is the same as that for a direct investment in the stock). Ignoring possible 
external constraints on time frames and margin requirements, this question is entirely a 
question of the risk profile the investor wishes to assume. As a first quantitative measure 
of the risk we would suggest to look at sd(P)/E[P], i.e. the standard deviation of the 
option value per unit price. In the first graph in Figure H] we plot this, for a stock with 
r = 0.02, cr = 0.25, Sq = 25, as a function of K and for a range of values of T. As we 
might expect, out-of-the- money options (K < Sq) are apparently much more risky than 
in-the-money options (K > Sq), and whereas for deep in-the-money options shorter term 
means lower risk, for deep out-of-the-money shorter term means greater risk. (The behavior 
close-to-the-money is interesting, as the insert in the figure illustrates.) 

From the exact formulae for the variance of the European put option it can be shown 
that as K tends to infinity the variance Var(P) tends to S^ie? T — 1), which is the variance 
of the stock price (discounted to current values). Thus the first plot in Figure [4] shows 
the evidently sensible result that for fixed T the risk associated with buying put options 
decreases with the strike K, ultimately decreasing to the risk of the stock itself. 

Another quantity the investor may wish to look at is the PEW. This is also displayed 
in Figure HJ In fact we see that PEW behaves very similarly to sd(P)/E[P]. Indeed, by 
Chebyshev's inequality 



In practice, as usual with Chebyshev's inequality, this is not a very useful inequality quan- 
titatively, but it does, however, correctly indicate that the general trends of PEW and 
sd(P)/E[P] are similar. 




and thus 



PEW 



P(P = 0) < P(|P-E[P]| > E[P]) < 
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Figure 4: sd(P)/E[P] and PEW for European put options on a stock with r = 0.02, 
o = 0.25, S = 25. Includes blow-ups of the intersection regions. 

It should be emphasized that the risk we are discussing here is risk as viewed from the 
buyer's side. For PEW this is clear; for the buyer of the option, high PEW means a high 
probability of losing 100% of the investment, which is an uncontroversial, if rather naive, 
measure of risk. For the seller the opposite is true; high PEW is a measure of safety for the 
seller. This reflects the general asymmetry between buyer and seller in matters of risk — a 
buyer prefers to buy low risk, a seller prefers to sell away high risk. For sd(P) /E[P] things are 
a little less clear. If we are discussing an options trader, who is buying (or selling) an option, 
say, to add to (or remove from) a portfolio, then indeed, sd(P)/E[P], which measures the 
spread of the return from the option, is a measure of buyer's risk, and the same asymmetry 
between buyer and seller exists. (Note, however, that if the trader holds other options on 
the same underlying then it is of course a mistake to consider just the risk on an individual 
asset in the portfolio and ignore correlations.) The potentially confusing situation comes in 
considering an investor who wishes to write, say, a single in-the-money put option, with a 
view to actually owning the stock at a later date when the option is assigned. Despite being 
a seller of the option, this investor may well prefer lower sd(P)/E[P]. For this seller, the 
sale does not eliminate risk, but adds to it — so the considerations of a buyer may be more 
appropriate. 

6 Risk-averse/Risk-seeking Option Pricing. 

The previous section concerns the use of risk measures of options in individual investment 
decisions. In this section we consider the following question: Is it possible that the market 
incorporates a measure of risk in option prices? The simplest imaginable way to incorporate 
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Figure 5: Effective volatility of European call and put options priced using (fT4"|) . as a function 
of K. Values of q from top to bottom: —0.02, —0.01, 0, 0.01, 0.02. Stock parameters r = 0.02, 
o = 0.25, So = 25; option expiration T = 0.5. 

risk premium is to consider a pricing formula of the form 

option price = expected value — q x standard deviation (14) 

where q is a parameter quantifying the average risk aversion of option buyers. Since in 
the classical Black-Scholes formulas option prices are monotonic increasing functions of 
volatility, we can translate the above change in option prices into effective volatility values. 
Note that for q > (net risk aversion) the above formula gives prices below classical Black- 
Scholes, and hence effective volatility is lower than input volatility. For q < (net risk 
seeking) the reverse is true. In Figure |5] we plot the effective volatility of European call 
and put options as a function of K for a few values of q. We assume an underlying with 
r = 0.02, a = 0.25, So = 25 and look at options with expiration T = 0.5. 

We see that for q < our model predicts a volatility smile and for q > a volatility frown. 
The plots in Figure |5] are somewhat misleading, as very large ranges of strikes are shown 
(differing in the 2 subplots). For deep out-of-the-money and deep-in-the-money options, 
the classical Black-Scholes price is very insenstive to the volatility. Thus a small change in 
price requires a substantial change in volatility. This is the reason for the sharp upturns and 
downturns in effective volatility, particularly noticeable in the case of the call. Nevertheless, 
it seems that Equation ffT4l provides a 1-parameter extension of the Black-Scholes model 
which can explain a volatility smile. The standard extensions — local volatility models 
[TJ [TO] and stochastic volatility models (THJ [UJ [13] — typically need at least 3 parameters. 

The term structure of effective volatility implied by Equation ( fT4"l) is complicated and 
we do not describe it here. However, we note that since q is a parameter quantifying risk 
aversion, it is quite reasonable to take q to depend on the timeframe T, thereby allowing 
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any desired term structure to be built into the model. As of yet we have not tested to see 
whether the model can provide a reasonable fit to data. 

7 Concluding remarks 

In this paper we have explained methods to calculate the variance (and the PEW) of stan- 
dard options, and outlined potential applications. Despite the fact that we have not yet 
performed comprehensive empirical studies, we think the potential usefulness of the variance 
and/or the PEW is self-evident. A further point may help to make this clear: In the classical 
Black-Scholes option formulas, the price of an option depends on 5 parameters So, K, T, r, a. 
By rescaling units of time and money 2 of these parameters can be eliminated. By moving 
to a frame in which all money accumulates at a rate r, dependence on r can be eliminated. 
Thus there is only non-trivial dependence on 2 effective parameters, which are typically 
taken to be the combinations ay/T and K f° r r ■ If an option is priced using the classical 
Black-Scholes formulas, then two pieces of data are being compressed into one — in other 
words, information is being lost. By also looking at the variance we take all information 
into account. 

In addition to empirical studies and investigation of potential applications, there are 
many extensions of the current work that come to mind, and we mention just a few: 

• The variance formulas for European options in this paper should be extended to cover 
the case of options on (continuous) dividend-paying stocks, for which exact formulas 
should still be available. 

• For American options, we have seen exact formulas are not available, but we believe 
it should be easy to provide good approximations, particularly in the limit of low 
volatility, following the work of Widdicks et al. [20J |8] 

• Similarly for options on multiple assets, asymptotic approximations should be avail- 
able. Good methods for options on multiple assets should involve incoporating asymp- 
totic results into numerical schemes to have the best of both worlds. 

Acknowledgements. We thank Yaniv Zaks and Evgenia Apartsin for useful conversations. 
This work includes material included in the M.Sc. thesis of the first author. 
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